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In this paper, I have used direct methods in order to evaluate the inflationary dynamics of f(R, OR) gravity 
with Lagrangian /(/?,□/?) = R + c R 2 + 2™i CjRO'R. This is found to have equivalent solutions to R + R 2 
gravity, and to be dynamically equivalent to chaotic inflation. Perturbation spectra are calculated, and fitted to 
the COBE normalisation. This is found to correspond to a large field inflation model. 
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I. INTRODUCTION 

Ever since Einstein and Hilbert first noticed that Einstein's 
Field equations could be derived from an action principle, 
there has been interest in the effects of inserting additional 
terms into the Einstein-Hilbert action. However, it was the 
rise of the inflationary paradigm [ 1] [2] [3], coupled with the 
discovery of field theory, that has provided motivation for a 
more widespread investigation of this field. The first and fore- 
most motivation is that we expect that gravity should have a 
field theoretic description, in which case the Einstein-Hilbert 
action should be treated as some kind of effective field theory. 
In this case we expect that there should be present a very gen- 
eral class of terms with higher powers of the curvature tensor, 
that are suppressed according to some mass scale. It is hoped 
that this will also make gravity into a renormalisable theory at 
low energy. These terms can have significant effects either at 
early times when the universe has higher curvature, or at very 
large scales where tiny effects can become cumulatively sig- 
nificant, at least in areas where the mass density is sufficiently 
low. It is generally thought that additions based around the 
Ricci Scalar should be effective toy models for more general 
tensor additions. 

Secondly, the question of finding a theoretically well moti- 
vated construction for inflation is still open. Within the slow 
roll approximation, scalar fields can be fitted to the data with 
an arbitrary degree of precision, fitting one into a well mo- 
tivated particle physics construction has proved to be much 
more difficult. Examples exist in the Minimal Supersymmet- 
ric Standard Model or the Higgs Inflation Model 

|0, but until experimental verification, they remain only theo- 
ries. Until such questions are answered, positing the existence 
of a slowly rolling scalar field seems no more likely (and in 
many ways less so) than the idea that gravity could itself drive 
inflation||8|]. 

The aim of this paper is to investigate the inflationary dy- 
namics of adding a series of arbitrarily high derivatives of the 
form RD k R. The addition of an infinite series of such terms 
was previously found to be ghost free, and to give rise to an 
asymptotically free theory of gravity, and a bouncing cosmol- 
ogy m- 



II. DERIVING THE EVOLUTION OF THE BACKGROUND 

In this paper we consider a Lagrangian given by 

L = /(/?, OR) V=S, CD 



where f(R) is given by 

oo 

f(R) =R + c R 2 + Y aROR, 



(2) 



and □ has its usual meaning as the d'Alembertian. We use 
here a result of Schmidt 1 10] for the Lagrangian ([TJ giving 

1 6L fe ij 
—— =GR' J - — - G' ,J + g J DG 



A=\ A=l 

(3) 

where we have defined the following functions: 



j=A 



l df 
dDJR' 



G = F . 



(4) 



If we assume that the dimension of our space time is 4, then 
we can take the trace of the variational derivative to find 



3DG = 2/ - Gfl - 2 ]T F A D A R + 2 ^ Fa^D^ 1 R)' J . (5) 



A=\ 



A=\ 



As is usual in a homogenous isotropic universe, only the 00 
component of the field equations and the trace provide inde- 
pendent constraints, and it will be sufficient to work with the 
much less complicated trace equation to establish the evolu- 
tion of the curvature scalar. 



A. Finding solutions 

Before inserting our Lagrangian into the trace equation, it 
is instructive to input two simpler Lagrangians defined by 

f(R, OR) =R + aR 2 , (6) 
f(R, OR) = R + aR 2 + ROR. (7) 

In the first case the trace equation is reduced to 

3DG = 2/ - GR, (8) 

with 

G = % = 1 + 2aR, (9) 
dR 
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and consequently we have the equation of motion for R as 



DR = —R. 

6a 



(10) 



Under the assumptions of homogeneity and isotropy the 
d' Alembertian reduces to a double time derivative (clearly un- 
der such conditions all the space derivatives must be zero), in 
which case it is clear that R has the solution 



R = A exp(z 



V6a 



(11) 



It is interesting that the sign of alpha gives the difference be- 
tween oscillatory and exponential behaviour. Notice also that 
R — is a solution, that is, a flat space time remains a stable 
solution in the absence of matter, and indeed, is an attractor of 
the vacuum solution. The equation of motion does not spec- 
ify whether the curvature is increasing or decreasing, we wish 
to study inflationary modes, and therefore choose that a > 
and R < 0. Using again the assumptions of homogeneity and 
isotropy we are left with a Robertson- Walker space time, and 
this implies that R = 12H 2 HI. With (TO]) given above, and 
the usual proscription for H in terms of the scale factor, it is 
clear that this relation can be integrated to yield that 



lna(f) = -A V2aexp( — =) + B, 



yielding 



a(f) = B' exp(-A V2aexp(-- 



=)), 



(12) 



(13) 



where B/B' have been introduced as constants of integration. 

Let us continue with the calculation of the background evo- 
lution of 

f(R, OR) = R + aR 2 + yRDR. (14) 

Using the trace of the field equation (0, we find that 

6aUR + 6yD 2 R = R- 2yRUR + 2yViRV% (15) 

and under the conditions of homogeneity and isotropy we can 
make the identification that V = j t , leading to 



2y 



2\ 



dR\ 

at) j 



One can immediately see that for any exponential solution the 
left hand side vanishes identically, so let us put R = expBt 
into the right hand side to produce, 



= 6yB A + 6aB 2 - 1, 
which clearly has solutions if 

-6a ± ^36a 2 + 2Ay 



Yly 



(17) 



(18) 



This clearly has real solutions for B if y > 0, or if y < 



but \y\ < \a 2 . Given that we expect each term with higher 



derivatives of the metric to be associated with higher powers 
of some mass scale M, combined with some coupling constant 
of order 1, we find that the second of these conditions is nearly 
automatically satisfied, as it gives 



I'M < 



(19) 



Thus we see that despite the extra terms the solutions are 
nearly identical, still yielding exponential behaviour of R, 
with the corresponding form for the scale factor Eq.([T3l). 
There is a pleasing symmetry about this, that extra terms do 
not disrupt the form of the equations. 



B. Solutions to the full infinite derivative Lagrangian 



In this section we present results based on a direct analysis 



of 



f(R) =R + c R 2 + Y aRHR, 



(20) 



and we need to insert this into the trace of the field equation 
Eq.©. This is a slightly tedious calculation made easier by 
noting 2F - GR = R is an easily established identity. In order 
to see the nature of the solutions it is convenient to group the 
terms according to whether they are linear or quadratic in R, 
giving 

3nG-R = Y(v i F A v i n A ~ 1 R-F A n A R). (21) 



It is easy to see that an exponential solution can only solve this 
if one of the sides of this equation reduces to zero. If we make 
the assumption of homogeneity, so that we can drop the space 
derivatives, then we find immediately that the right hand side 
will reduce to zero for any exponential solution. In order to 
understand the solutions physically, we shall make explicit the 
mass dependence at this time, by writing, 



R = A exp(-fcMf) 



and, 



ypi+2 ■ 



Using these, we are left with 

= 3D 
reducing to 



i +2 ^ +2 y 1 

M 2 4-? 

1=1 



(22) 



(23) 



(24) 







6^ib 2i+2 -h 

i=0 



(25) 



upon insertion of the trial solution. This type of polynomial 
will have (at least one) real solution for b under some very 
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general constraints. In particular if we consider that all the 
Aj are about 1 , and then this is a becomes a simple geometric 
series in terms of b. In particular we can use the geometric 
sum noting that, 



1 -b 2 ' 



(26) 



giving b = -4= . We had at this point a free choice of sign under 
the square root, but we elected for the positive sign given that 
we are looking for inflationary solutions. At any rate, negative 
solutions will be unstable as the curvature will tend exponen- 
tially towards a singularity. In practice the exact value of b is 
not necessary, given that we shall treat the mass scale M as a 
free parameter to be determined by the COBE normalisation. 



III. INFLATIONARY DYNAMICS. 



the extra terms of the Lagrangian all reduce to the form of an 
R 2 term during the inflationary scenario, when homogeneity 
and isotropy can be invoked, this indicates that our theory has 
a single extra (scalar) degree of freedom compared to stan- 
dard Einstein gravity. It is then no surprise to find a dynami- 
cal equivalence to Einstein Gravity plus one scalar field. This 
equivalence has also been extensively investigated in the form 
of a conformal transform 1 . In order to show the equivalence of 
our theory with scalar field theories, we shall use the slow roll 
approximations to define an implicit field <p, which we shall 
then show fulfils the canonical scalar field equation, justify- 
ing the use of the slow roll predictions. As in the use of the 
conformal transformation, the field <p is a purely mathemati- 
cal convenience - in practice it is a substitution for R, reducing 
our field equations to the form of a scalar field. However, un- 
like a conformal transform, there is no mixing of degrees of 
freedom. 

In essence we use the relations, 



In this section we shall investigate the inflationary dynam- 
ics, in particular we shall find an exact form for the scale fac- 
tor, and show that inflation will end naturally. We will then 
show that the the theory is slow rolling, and that it shares a 
dynamical equivalence with a scalar field theory, which shall 
be used to evaluate the spectrum of density perturbations in 
the following section. 

Firstly, we note that under the assumptions of homogeneity 
and isotropy, the Freidmann-Robertson- Walker metric is nec- 
essarily a solution to any metric theory of gravity, and thus 
we can make use of the well known relation (see for example 
lfl2tl ) R = 12H 2 . It is not difficult to integrate this to find that 



(27) 



3M 2 PI H 2 = V(cf>(t)) 



a(t) = «o exp 




The subscript indicates that the quantity is defined at t = 0, 
which does not necessarily correspond to the beginning of the 
universe, but to an arbitrary time when we choose to define 
both R{) and oq. We shall use that Rq = M 2 pl as this seems 
to be a natural scale for R at the start of the inflationary era. 
We also note that despite the odd form of a, it is a quasi de- 
Sitter space, at least initially, as provided -4p << 1, then 
exp(^p^) a 1 - ^y±. With this exact solution known, we can 
investigate the behaviour of a( t), and we find that 



a(t) — a -i/ — 




bM) I bMt^ 

— j ex T— 

(28) 

is monotonically decreasing, and reaches a — 0, correspond- 
ing to the end of inflation, at a finite time given by 



3//0 = - 



dVW)) 
d<f> 



leading to 




l6 2 M 2 pl 



(30) 
(31) 
(32) 



(33) 
(34) 



This method of direct equivalence should be accurate, pro- 
vided that (f> fulfils the canonical wave equation to within a 
good approximation, that is, ip « 3H^>. It is easy to check 
that this is certainly satisfied provided that <p is greater than 
around 5Mp/. This is important, as it ensures that the pertur- 
bations will have the same classical evolution as in the slow 
roll paradigm, when cosmologically significant scales move 
outside the horizon. It is generally believed that the effects of 
a truly quantised gravity will be sub-dominant to the effects of 
the classical evolution. Thus this extended gravity theory will 
have cosmological predictions that are necessarily similar to 
the slow roll paradigm. We also note that the form of the po- 
tential is essentially that of chaotic inflation, but the parameter 
A is now seen to be a function of the mass scale M, potentially 
explaining the otherwise unnatural smallness of the coupling 
constant. 



tend - ~T~T7 m 

bM 




(29) 



If we look at the parameter jj? then we find that it is less than 
1 until exactly t em {. Physically, this parameter is equivalent to 
the slow roll parameter of a scalar field theory, (see, for exam- 
ple, lfl3ll ). This should come as no surprise, as in this theory 



1 I believe that conformal transformations may obscure physically relevant 
similarities between theories. In particular it is not always clear what ap- 
proximations made in the conformal frame correspond to, when a con- 
formal transformation represents a point-wise scaling of such physically 
meaningful quantities as mass-energy and length. The mathematical equiv- 
alence of two frames related by a conformal transformation is only pre- 
served under the redefinition of units in the conformal frame to account for 
the point wise rescaling of length, see for example, fl4ll 
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IV. COMPUTATION OF THE SLOW ROLL PREDICTIONS 



Firstly we begin with the slow roll parameters, 

£ ■ 



>1 



— 

T 



8M 2 pl 



12M% 



(35) 
(36) 



Thus we can see that inflation comes to an end at (f> — 
2 y2Mpi, which is precisely t en d, given above. The next im- 
portant result is the total number of e-folds, this is given by 



N = In 



a(t m „) 



(37) 



= (m VS (exV^stan ~ exp(-f^ t end j) . (38) 

= li(^-€j ( 39 ) 

Thus the dominating factor in the number of e-folds is 
where inflation starts, rather than when it ends, although this 
is not altogether surprising given that inflation effectively ends 
at a fixed time controlled only by the parameter M. In prac- 
tice the choice of t stan is implicitly tied to the choice of Rq, it 
is convenient therefore to choose t — as the start time. Given 
that it is expected that Rq will be of the order M 2 pl and that M 
will be several orders of magnitude smaller, it can be seen that 
this theory produces more than enough inflation to account 
for observation. We will check the value of M against the 
COBE normalisation next. We use the simplifying assump- 
tion of [13] that cosmologically interesting scales cross the 
horizon when N rema i n i ng = 50, 2 corresponding to <p = 2Q.2M P i. 
We continue the paradigm by evaluating the spectral index n 
and its running, evaluated at Nsq. (All data for reference has 
been taken from lfl5ll ). 



n = 1 + -6e + 2T] = 0.941 
-f&r = 5.77 x 1(T 4 . 



(40) 
(41) 



We are now in a position to evaluate the COBE normalisation, 
which is given by the expression 

1 V 



<w = 



\5Qn 2 M A pl e ' 



(42) 



As usual this is evaluated at the pivot scale, corresponding 
to around 50-e folds before the end of inflation, giving the 
expression 

b 2 M 2 



5 l H {k) = 67.2 x 



M 2 



(43) 



which, using b — -^= as calculated in section 1, and the stan- 
dard COBE normalisation we find that M = 6.16 x 10" 6 M W = 



1.5 x 10 l3 GeV/c 2 . This is in keeping with results obtained via 
the conformal transformation, see for example, 11611 . The final 
piece of analysis concerns the scalar to tensor ratio, given by 



r = UAe, = 0.24. 



(44) 



In reality, this is a little high, the WMAP data strongly favour- 
ing r < 0.2. However, as alluded to in an earlier foot note, this 
is not actually an insurmountable obstacle, as the choice of 
Nremaining is based on the thermal history, with many authors 
preferring 55 rather than 50, and some textbooks suggest that 
values as high as 70 are possible, if disfavoured. It is per- 
haps sensible to think that ( fT3l may govern the evolution of H 
for around one Hubble time before the radiation density has 
grown enough to dominate, this would tend to increase N* by 
a few e-folds. However, a choice of 70 is needed to bring the 
results into the 2-cr region of the WMAP probe, and that is 
probably too high, unless a theory of reheating should cast the 
standard thermal history of the universe in doubt. 

We have not discussed reheating at all in this model. The 
scalar field above does not correspond to a real scalar field, 
and will not oscillate around a minimum as in the standard 
scalar field examples. However, at the point where accelerated 
expansion ends, there is still a large energy density locked in 
the Ricci scalar, and if this couples to some particle physics 
model, as in the Higgs [7] inflation model, then this could de- 
cay into the particle sector. Alternatively, the particle physics 
could be wholly decoupled from the gravity action, and con- 
tain fields of its own which can decay into the standard model. 



V. CONCLUSION 

We have shown, via a direct evaluation, that our infinite 
derivative Lagrangian has the same inflationary solutions as is 
found in R + R 2 gravity, and that these solutions are function- 
ally identical to that of a scalar field with a 0-fourth potential 
when cosmologically interesting scales pass outside the hori- 
zon. This author eschewed the use of a conformal transforma- 
tion in this example, seeing it as an unnecessary complication, 
and notes that here the use of a direct method has revealed the 
phenomenological equivalence with a phi fourth scalar field 
which is obscured (though presumably it still exists - see 11411 ) 
by a conformal transformation. This therefore has similar in- 
flationary dynamics to chaotic inflation, only now the small 
value of the dimensionless coupling constant lambda (identi- 
fied with the constants in Eq.d34"l)) needed to match the COBE 
normalisation is seen to be motivated as it is found to be de- 
pendent on the mass scale M. Unfortunately, it therefore suf- 
fers the same fate as phi-fourth inflation, which is generally 
thought to be outside the 2 sigma parameter space of the 5 
year WMAP data. 
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2 This value corresponds to the generally accepted number for when observ- 
able scales passed outside the Horizon, but it depends on the assumption 
of Einstein gravity. Even then it is only an approximation and so it is pos- 
sible to increase this number as high as 60 without any serious constraints. 
Further, given that the behaviour of H becomes non canonical a few e-folds 



before the end of inflation, and until reheating produces sufficient radiation 
to enter radiation domination, the faster than expected (cf radiation domi- 
nation) decrease in H should could leave us with a higher than expected N, 
and could favour N=60. 
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